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Abstract 

We discuss some aspects of the search for CP asymmetry in the three body B decays, revealed 
through the interference among neighbor resonances in the Dahtz plot. We propose a competitive 
method to extract the CKM 7 angle combining Dalitz plot amplitude analysis of — > K^TT'^Tr~ 
and untagged B^, B^ Kgir^TT^ . The method also obtains the ratio and phase difference between 
the tree and penguin contributions from B^ and B^ K*^tt^ decays and the CP asymmetry 
between B^ and B^. From Monte Carlo studies of lOOK events for the neutral mesons, we show 
the possibility of measuring 7. 

PACS numbers: 
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I. INTRODUCTION 



According to the Standard Model, CP violation is only possible through the CKM com- 
plex parameters. Interference between processes with different weak phases contributing to 
a same final state can generate an asymmetry in charge conjugate B meson decays. The 
most established methods to extract the CKM phases exploits the interference between the 
phases generated in the B^ oscillation and the phase of the decay of these particles in 
the same final state. Using such approaches on the B^ J/'^Ks decay, the CKM f3 angle is 
measured to a good precision. However, the other two CKM angles, a and 7, retain sizable 
uncertainties, due to both experimental and theoretical difficulties. LHC should determine 
7 to ~ 5° in one year of data taking[l|. 

Any asymmetry observed in charge conjugated B decays is interpreted as a manifestation 
of the CKM phases. In the two body decay scenario, direct CP asymmetry {Acp ) was 
established by the experiments BaBar [2] and Belle [sl by simply counting the difference in 
number of events of the 5° — K^tt~ and B^ —>■ K~7t^ decays. In three body B meson 
decays, BaBar and Belle claim the observation of direct Acp in the channel B 

X±pO measured through the amplitude analysis of the final state B^ K^ti^h^ . They see 
a difference in the decay fractions among the intermediate states in the B~^ and B^ samples. 
That is the difference between the square of the fitted amplitude for the intermediate state 
integrated on the phase space. This quantity is a straightforward application of the two body 
strategy into the three body decays, including the requirement of a strong phase difference 
to observe Acp ■ 

Beyond the fraction asymmetry quoted by Belle and BaBar, one could explore the asym- 
metry associated with the interference among two intermediate states, i.e. the interference 
between two neighbor resonances decaying to the same three body final state. As will be 
discussed in detail in section III, whenever the phase difference between two resonant ampli- 
tudes involves the weak phase, one should expect to see Acp in the interference terms of the 
Dalitz plot distribution. This contribution to the asymmetry depends only on the existence 
of a weak phase difference and it can be observed even in situations where there is no strong 
phase difference. There are a lot of opportunities for investigation of the asymmetry in the 



Dalitz plot distribution rather than difference in the decay fractions alone 



a. 



The idea of extracting 7 exploring the Dalitz interference among intermediate amplitudes 
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in three body B decays was initially proposed for 71^71"*" vr" [s^, where Xco plays a 

fundamental role as reference channel. However, the method is statistically limited by the 
low contribution of the Cabibbo suppressed amplitude B^ 



XcoTi"^ isl, and also estimated 

through the already observed Cabibbo allowed B^ —>■ XcoK^ channel {4, 5]. Recently, a 
new effort has been made to study 7 in S ^ K-k-k decays, however this approach requires 
time-dependent analysis 10|, having to deal with tagging inefficiencies that could be as 
large as 90% at LHCb [jj. 

We present a new possibility to measure the CKM 7 angle, through the interference terms 
in the amplitude analysis of three body B decays. The method combines the experimental 
data from B^ K^-n-^n^ and B^, B^ KsTt^tt^ . From the charged B decays we have the 
penguin amplitude coefficient and phase of the channel K*tt. We then use this result in an 
untagged analysis of the neutral system, to extract the tree component strong phase and 
amplitude coefficient and the weak phase 7. The method is based on the ability to measure 
independently the relative amplitudes and phases for B^ and B^ decays in a joint untagged 
sample. 



II. AMPLITUDE ANALYSIS 



The standard amplitude analysis, extensively used in non-leptonic three body heavy 
quark meson decays, is based on the Isobaric model. In these analyses, the two body 
resonant intermediate state amplitudes are represented by Breit-Wigner functions multiplied 
by angular distributions associated with the spin of the resonance. The various contributions 
are combined in a coherent sum with complex coefficients that are extracted from maximum 
likelihood fits to the data. The absolute value of the coefficients, normally refereed to as 
magnitudes, are related to the relative fraction of each intermediate state and the phase takes 
into account rescattering or the final state interaction (FSI) between the resonance and the 
bachelor particle. This phase is considered constant because it depends only on the total 
energy of the system, i.e. the heavy meson mass. The presence of strong final state constant 
phases can be directly revealed through the interference between two different amplitudes. 
Thus, the total Isobaric amplitude, written as functions of the Dalitz plot variables, S12 and 
Si3 is: 
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N 

At = a^Re^^^«^7Vi?(si2, sis) + J2 ^ne'^' -^nisu, S13) • (1) 

n=l 

The fit parameters are the coefficient magnitudes, a, and the phases, S. The non-resonant 
amplitude, Anr is usually represented by a constant. Each resonant amplitude, An{n > 1), 
is written as a product of three terms: 

An= 'FnX'MnX BWn (2) 



The first term are the Blatt-Weisskopf damping factors 13|. "^Ain is a term which accounts 
for angular-momentum conservation and depends on the spin, J, of n^^ the resonance, such 
that it is 1 for a spin-0, — 2 | ps || p2 | cos9 for spin-1, and |(| pa || p2 \y{3cos^9 — 1) for 
spin-2, in the case of the resonance formed by the pair of particles 1 and 3. p2 and pa are 
the the 3-momentum of particles 2 and 3, and 6 is the angle between particles 2 and 3, all 
measured in the resonance rest frame. The last term is a relativistic Breit-Wigner function 
given by: 

BW^ = ^ . y (3) 

where m„ is the resonance mass and r„(sjj) is the mass dependent width. 

The relative contribution of each intermediate state is given by fractions defined as: 



fn - r.' :\:'\ :l • (4) 





^nAn 
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where / dsi2dsi3 J2jk I ^i-^i'^k-^l I related to the total number of events of the non-leptonic 
three body heavy meson decay. 



III. THE DIRECT DALITZ PLOT ASYMMETRY -Adp 

One important feature of the Dalitz plot associated with the interference between two 
neighboring resonances is that there will be an Acp manifestation, independent of the strong 
phase differences. To illustrate this point, let us develop a simple example of two resonant 
scalar amplitudes, interfering in a particular region at the three body decay phase space. 
Suppose one of the resonant amplitudes, Ai, has a tree component, with magnitude and 
strong phase 5t and weak phase 7, plus a penguin contribution with ap and 5p and no weak 
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phase, interfering with another resonant amphtude, Ar, with magnitude and phase and 
6r and the absence of weak phase. The charge conjugate total decay amphtudes are: 



A- = [are''^^^""'^ + ape'^'']Ai + aRe'^'^AR 



(5) 
(6) 



and 



|^+|2 _ 1^ |2 ^ -4aTapsin7sin(5r - ^p)|AP 
-AaRttT sin 7(sin {St - Sr) UeiAiAl) + cos {6t - Sr) ImiAiAD) (7) 

From Equation 7, we see that the first term has a similar behavior as the two body decays, 
i.e. Acp depends on the magnitude of the tree component and also on the phase difference 

— 5p), which must be non zero for Acp to be observable. This is the term responsible for 
;he Acp quoted through the difference in integrated fractions made by Babar and Belle ^ 
^. However the presence of second or third terms, the interference ones, depends only on 
the magnitude of the tree component, since the phase difference has a sine form associated 
with the real part of the BW's product, but also a cosine component associated with the 
imaginary part of the BW's product. So the interference will always produce Acp provided 
that the magnitude associated with the weak phase is significant. 

Interferences in the Dalitz plot do not allow to associate directly the fraction, defined 
in Eq. 4, with the number of events in a intermediate state channel. The relative phase 
between two amplitudes sharing the same phase space region can produce destructive or 
constructive interferences, decreasing or increasing the number of events in the three body 
final sample, in which case we lose deterministic information of the source of the events in 
the region. Due this difficulty to the concept of the direct Acp to the three body decays, 
the asymmetry observed from fractions would be better refered as quasi- direct Acp ■ 

The direct Acp in the three body decays, in the sense of event counting, would be 
more conveniently defined by the notion of two asymmetries, one local and one global. A 
global Acp is the difference between the total number of events of two conjugate three body 
decays. The local would be the difference in the number of events depending on the Dalitz 
variables Sij. In general, the global asymmetry loses part of dynamics information, since it 
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is the integral over the phase space. On the other hand, local asymmetry accounts for the 
difference of the magnitude square, or the first term of Equation 7, and also the asymmetry 
coming from the interference terms among two resonances, or the last terms of Equation 7. 

The local CP asymmetry of the charged B meson three body decay can be defined through 
the subtraction of the two Dalitz surfaces relative to the charged conjugated B signals: 

. ^ A^^(-Sl2,-Sl3) - N-{si2,Si3) 

N+{s,2,s,,) + N-{su,s,s) ^ ' 

we call Adp the local asymmetry at the Dalitz plot, with N^{si2, S13) being the number of 

events in S12, S13. 

For the neutral systems with and B'^ decaying to the same final state, Kstt^tt^, Kspp, 
7r°7r+7r^ for example, Burdman and Donoghue 6| propose a convenient definition of local 
Adp- They showed that for a decay with \p/q\ = 1, the asymmetry is time independent and 
can be written as: 

, _ iV°(gl2,gl3)-iV°(gi3,gl2) 

A measure of local Adp is, obtained by comparing the number of events in points of the 
Dalitz plot symmetrical with respect to the diagonal line S12 = S13. 

To quantify the total amount of asymmetry present in a three body decay, one can 
compute the sum of the modulus of Adp, taking in account properly statistical fluctuations. 

IV. METHOD 

We now introduce a new possibility to measure 7, based on a time-independent analysis of 
B Ktiti decays. The most relevant observed intermediate states for B^ ii'^7r+7r~and 

are summarised in Tabled with the respective tree and penguin 
contributions, and CKM weak phases contained in the amplitude. 

Based on SU(2) flavor symmetry, we expect the same penguin amplitudes for the four 
processes 5^ K*ti^ and B^,B^ — >■ K*'^7T^. This property allows to extract the B"^ 



K*TT^ penguin parameters 231], and use them in the B^ and B^ joint fit to measure the tree 
phases, from which 7 is obtained. 

Since we measure only relative magnitudes and phases we must elect one resonance to 
have fixed parameters. The compatibility between penguin parameters from B^ and B^, B^ 
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TABLE I: Resonances and dominant contributions from B Ktttt decays . We denote allowed, 
suppressed by color and penguin amplitudes by T*^, Ts and P. 

to i^r*7r, requires that the amphtude analysis be made relative to an anchor resonance that 
has the same amplitude for B charged and neutral. For that we chose the non CP violating 



B ^ K amplitude. The asymmetry measured by Belle [15| for the channel B"^ K'^Xci 
is Acp = —0.01 ± 0.03 ± 0.02, indicating that the dominant contribution is a tree diagram 
without weak phase. 

Following is a schematic representation of the method, with the arrows pointing its flow: 



B° ^ K*+n- : ap e'^^ + e^^- 
B^^K*-n+ : apc'^^ + aj.e'^~ 

7 = ■ (10) 

where 9"^ = {St ± 7). We underlined the parameters measured in each amplitude analysis. 



The method is based on three basic and well accepted hypothesis that can be tested 
during the analysis. Hypothesis 1: the dominant contribution of B"^ K*t[^ is VbtV*^ P. 
This has been indicated by the BaBar result [4]. The test is to check if the K*tt^ 
and B^ K*^^ amplitudes are the same. Hypothesis 2: the penguin components from 
B^ K*7T^ and B^,B^ K*^-!!^ are equal [igI . \y\ . Hypothesis 3: Xco have the same 
amplitude for B"^ — > K^Xco and B^,B^ KsXco- The experimental test for the second 
and third hypothesis, consists in the equality of the tree magnitudes measured from the 
intermediate process B^,B^ — * K*^7r^. The confidence in the result of 7 is based in the 
fulfillment of all hypotheses, whereas the failure of either implies interesting unexpected 
effects, like a big W annhilation contribution or even a new contribution for the B"^ — > K'^Xco 
other than the tree component. 



V. K^TT+TT- AND 5°, fiO ^ KsTT+n- AMPLITUDES 

To obtain the B^ K^n'^ii' parameters, we apply the usual maximum likeli- 
hood fit to the isobaric amplitudes given by Equation 2 with the intermediate states 
i = K*{890)°, K*{imf, p(770)°, /o(980) and K+Xco- 

Regarding the neutral system, separation between 5° and B^ to Kstt^tt^ is not trivial. 
The samples could be distinguished at the production level only through partner b particle 
tagging techniques. Furthermore, there is mixing among B^ and B^, introducing time 
dependence in the decay probabilities. 

The probability distribution for the final state Kstt^tt~ , independent of its origin, is given 
by the incoherent sum: |M(At)p + |M(At)p. The matrix element for observing the original 
B^ and B^ decay to the common final state are respectively: 



M{At) = e-(r/2-iA/)At cos {Am At/2) - iq/p A sin (Am At/2)], (11) 
M{At) = e-(r/2-^A/)At cos (Am At/2) - ^p/g^ sin (Am At/2)], (12) 

where A e A are time-independent decay amplitudes equivalent to Equation 1, but for 
B^ — > i^'sTT+TT" and B^ Kst^^tt' decays. 

Assuming the same production rate for B^ and B^, and in the case of \p/q\ = 1, it was 
shown [g], [1^ that this sum displays the interesting property of canceling mixing dependence 
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terms, 

|M(At)|2 + |M(At)|2 = e-^* ( + )• (13) 
Similarly to Equation 1, and are written as: 

A" = a^e' ^>^A^ + Y^ ate' At (14) 

i 

A^ = a^e' ^'^ A^ + Y^ cue' At (15) 

i 

As a fundamental step to achieve a 7 measurement, we propose a method of joint fitting 
the set of untagged plus B'^ time integrated event sample by doing a maximum likelihood 
fit to the Probability Distribution Function (PDF) given by 

I A0\2 I I J0|2 

PDF = L ' J , (16) 

where = J dsij dsjk, and iV° = J rfsjj dsjfc and a, 6 and a, 5 are free fit param- 
eters. 

In the standard single sample Dalitz plot amplitude analysis, described in section II, 
two parameters, one magnitude and one phase, should be fixed. What is measured are the 
relative contributions and phase differances. The magnitude sets the overall scale. This 
would be the case for each of the two charged modes with a total of four fixed parameters. 
In the joint fit, we have to fix three parameters, one magnitude, that sets the overall scale 
for both modes, and two phases, one for each c.c mode. In the joint fit, the freedom of all 
magnitudes of one of the modes, say B^, is what guaranties that we are able to have different 
total number of events for each mode. In our particular case we fixed a^, 6^ and 5^. We are 
able to measure a difference in number of events from and B^, searching for direct global 
asymmetric. With one single procedure it is possible to extract in an independent way the 
parameters from the amplitudes A^ and A^ and the ratio in number of events, calculated 
by the ratio of the normalizations N'^/N^. 

In the convention (^^5, 7r+, vr^) (1,2,3) for the final state particle numbering, the 
charge conjugation operation is equivalent to switching the Dalitz variables su ^ S13. The 
resonances K* bands are aligned in different axis for A^ and A'^, establishing some sort of 
signature for the event origin B^ or B^, as can be seen schematically in Fig. [H Although 
the TTTT resonances for 5° and B'^ overlap, the interference regions between Kn and vrvr 
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resonances are well separated. These interfence regions are vital for establishing the phases, 
and they guarantee that the fitting procedure finds a unique solution. Although we cannot 
distinguish events, the joint fit can identify two different overlapping surfaces relative to 5° 
and fio. 




FIG. 1: Schematic representation for the Dalitz plot of a)B^ i^svr+vr^, b)B^ — > ET^Tr+vr^ and c) 
the joint figure obtained for the untagged sample. Bands pictorially represent resonances regions. 
The vertical and horizontal bands are Kir resonances the in the diagonal are vrvr overlapping 
resonances. For simplicity we ignore angular effects in the distribution. Circles are displaying 
regions of interference within each or B'^ decay. There is no interfence between the horizontal 
and the vertical bands for they do not belong to the same mode. It is the non-overlapping of circles 
that assures the separation between both amplitudes in the joint fit analysis. 



VI. FEASIBILITY STUDY 

To investigate the feasibility and the error size in the joint analysis, we generated and 
fitted one hundred Monte Carlo samples of lOOK plus events, that can be obtained 
by LHCb experiment. We do not consider any theoretical or experimental systematic un- 
certainties. Detector and background effects are as well not taken into account. 

The parameters we used to generate the Monte Carlo samples are listed in Table [TTl 
column input. They are close to the values obtained in measured parameters for the B^ 
K^n'^Ti" decays . For the B^ and B^ — ^ Kgn^n decays we used the same parameters 
as in the charged mode, except for including a tree contribution, discussed below, for A'* (890) 
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and K^{14:30)'^ and the 7 phase, taken as 69°. The Dahtz plot distribution for one generated 
experiment of lOOK events is displayed in Fig. [21 where the small contribution of the 
resonance Xco can be seen. 




m 



Ks II- 



(GeV^) 



FIG. 2: B^, KsTT^TT events distribution from fast MC generated with table HIl parameters. 

One important issue concerning the 7 extraction, is the size of the ratio r = ax/cLp and 
the phase difference 6 = 6t — Sp from the tree and penguin amplitudes of the K* resonances. 
The smallest the value of r the more difficult it is to apply this procedure to measure 7. The 
theoretical knowledge of these quantities is model dependent. Some groups use factorization 
approach 17|], and obtain large r and small 6. On the other hand, non-factorisable approach 
for pseudoscalar-pseudoscalar B decay [19], presents an opposite scenario, with small r and 
large 6. The joint fit, applied to the real data, will be able to define which theoretical 
approach is more adequate, since it is possible to measure the parameters under discussion. 
For our study, we chose r = 0.45. 

As stated in section IV, the first step of the method would be to solve the charged mode 
Dalitz Plot, from which hypothesis 1 could be tested. The parameters ap and 6p would then 
be input to the neutral system fit. Since charged mode analysis presents no novelty and the 
statistical errors are expected to be low compared to the neutral decays we skip this step 
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KsTT+TT-) = 0.84 ± 0.12 



TABLE II: Monte Carlo simulation for B^, B^ KsTt^tt~ decay. We generated sample with the 
the parameters Oj and 6i for B^ and Oi and di for B^. The third column display the fit results with 
one hundred samples for each lOOK events. The last line is the ratio between B^,B^ KsTt^tt~ 
number of events. 

for this study. 

Each lOOK event sample of plus generated with the amplitudes and phases listed 
in the column input of Table [TTl were fitted with an unconstrained model with just three 
parameters fixed. In the last column of Table [ITl we display the central values of the Gaussian 
distributions of the 100 experiment's fit result and the errors are the width of the same 
Gaussian distribution. The total number of events of each B^ and B^ is the result of the 
integration of either amplitude within the Dalitz plot. The fact that the Xco magnitude for B^ 
is measured to be equal to the fixed ao is the test for hypothesis 3. The extracted quantities 
are in agreement with the generated ones and with small errors, stating the feasibility of the 
joint fit method. Using the parameters from Table HTl and the relations developed in section 
IV, we measure 7 = 69° ± 7°. 

If in the unconstrained analysis, the three hypotheses prove to be correct, we can reduce 
the numbers of free parameters assuming that the K* resonances parameters from B^ and 
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are given by the following tree and penguin sum: 

ae'^ = ap e'^^ + ar e'^^+^ , 

ae'~^ = ape'^^ +aTe'^^-\ (17) 

Instead of fitting for a,S,a and S, we fit for ap^Sx and 7 since we can fix ap,Sp from 
K^TT~^7T~ result. This way, constraining the fit with the hypothesis previously tested 
we reduce the error in 7. We measure with the constrained hypothesis, 7 = 69° ± 5° and 
also we obtained for the ratio between the tree and penguin component r = 0.45 ± 0.05 . 



VII. CONCLUSION 



We discussed some aspects of the search of Acp in three body B decays, emphasising the 
main characteristic of these decays, which is the interference among two neighbor resonances 
in the Dalitz plot. In particular, we discuss how the CKM 7 phase can be extracted applying 
the Isobar model 20|], extensively used in many charm meson three body decay amplitude 
analysis. We showed that the phase differences, coming from the amplitude analysis, can 
provide richer information than treating the three body decays as a straightforward extension 
of two body hadronic decays. We showed that if there is a weak phase difference between 
two intermediate state, Acp will be present in the Dalitz plot regardless of the strong phase 
difference. 

We presented a method to extract the CKM 7 angle using a combined Dalitz plot analysis 
from B"^ K'^'K^'K^ and 5°, B^ Kgn^TT^ . This method use three basic hypothesis that 
can be tested before one proceeds to the constrained fit. For measuring the B neutral 
analysis, we use a new technique of joint fit that allows us to extract, in an independent 
way, the amplitudes from two summed surfaces in a joint sample of B^ + B^ untagged 
events. We carried a simplified fast MC study to estimate the errors associated to the joint 
fit technique. Assuming a relatively low statistics for the anchor resonance Xco we obtained 
7 with a 5° error. 

During the analysis we can measure CP violation by counting the number of events from 
B^ and B^, information extracted in the joint fit, or exploring Dalitz symmetries as discussed 
in . We can also measure the ratio and phase difference from tree and penguin amplitudes 
of the K* resonance, defining which theoretical approach, factorisable or non-factorisable, 
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is more adequate. 

The methodpresented here is competitive with the other approaches to determine the 
CKM 7 angle [l| and needs the high statistics expected for the LHCb experiment, due to 
the small contribution of the reference channel B Y^K. However, in the case that the 
/o(980) resonance is dominated by the ss component 2l|], or even if the ratio between the 



tree and penguin is neghgible [22], the B — * fo{980)K amplitude could take the place of the 
charmonium reference channel in the analysis. 
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